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Abstract: TeV gamma-rays can result from the photo-de-excitation of PeV cosmic ray nuclei after their 
parents have undergone photo-disintegration in an environment of ultraviolet photons. This process is 
proposed as a candidate explanation of the recently discovered HESS source at the edge of Westerlund 2. 
The UV background is provided by Lyman-alpha emission within the rich O and B stellar environment. 
The HESS flux results if there is efficient acceleration at the source of lower energy nuclei. The require- 
ment that the Lorentz-boosted ultraviolet photons reach the Giant Dipole resonant energy (~ 20 MeV) 
implies a strong suppression of the gamma-ray spectrum compared to an E~ 2 behavior at energies be- 
low about 1 TeV. This suppression is not apparent in the lowest-energy Westerlund 2 datum, but will be 
probed by the upcoming GLAST mission. 

Two well-known mechanisms for generating TeV 
7-rays in astrophysical sources are the purely 
electromagnetic (EM) synchrotron emission and 
inverse Compton scattering, and the hadronic 
(PION) one in which 7-rays originate from ir° pro- 
duction and decay. Very recently, we highlighted a 
third dynamic which leads to TeV 7-rays: photo- 
disintegration of high-energy nuclei, followed by 
immediate photo-emission from the excited daugh- 
ter nuclei [1]. For brevity, we label the photo- 
nuclear process A + 7 — > A'* + X, followed by 
A'* — ► A' + 7-ray as "A*". Such a process may 
be operative in massive star formation regions with 
hot starlight. In this work we examine whether the 
A* -process could be the origin of the very ener- 
getic 7-rays recently observed, with the High En- 
ergy Stereoscopic System (HESS) of Atmospheric 
Cherenkov Telescopes, from the young stellar clus- 
ter Westerlund 2 [2] . 

RCW 49 is a luminous cloud of ionized hydrogen 
located towards the outer edge of the Carina arm, 
at a distance d « 8 kpc [3]. Embedded in RCW 49 



is the massive star formation region Westerlund 2, 
hosting an extraordinary ensemble of hot OB stars; 
presumably at least a dozen early-type O stars, 
100 B stars, and the remarkable Wolf-Rayet binary 
WR 20a [4]. For such a distance, the cluster core ~ 
5' results in a physical extent R ~ 6 pc. The total 
mass of the stars within this region is found to be 
rj 4500 M Q . The total wind luminosity of all these 
O type stars has been estimated as ~ 5 x 10 37 erg 
s~ 4 , and the stellar luminosity of known massive 
stars is L — 2.15 x 10 40 erg s~ 4 [5]. However, 
radio emission from the prominent giant HII re- 
gion RCW 49 requires a larger luminosity of ion- 
izing UV photons [6], Indeed observations us- 
ing the Infrared Array Camera (IRAC) on board 
the Spitzer Space Telescope indicate that the to- 
tal number of young stellar objects in this region 
is about 7000 [7]. Therefore, the above estimate 
of the total stellar luminosity should be taken as a 
lower bound. 

In the vicinity of WR20a, a clear excess of very 
high energy 7-rays was recently reported by the 
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H.E.S.S. Collaboration [2]. The significance of the 
excess is about 9a. Compared to the point spread 
function of the instrument, the source (termed 
HESS J1023-2013575) appears slightly extended, 
corresponding to an intrinsic size of the 7-ray 
source of about 0.2°; its center is slightly shifted 
compared to WR 20a. As expected for an extended 
source, the 7-ray flux is steady over time. 

By repeating the discussion of Cygnus OB2 [8], in 
what follows we obtain the expected 7-ray produc- 
tion through the A* -process in Westerlund 2. To 
compute the photo-disintegration rate of a highly 
relativistic nucleus (with energy E = AE^ = 
7 ArriN , where 7 is the Lorentz factor) on starlight 
per nucleon [9], 
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we must estimate the ambient photon distribution 
with energy spectrum n(e). In Eq. (1) we have ap- 
proximated the Giant Dipole resonant cross section 
by the single pole of the Narrow-Width Approxi- 
mation, 
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Here, e' is the photon energy in the rest frame of the 
nucleus, a /A = 1.45 x lCT 27 cm 2 , T = 8 MeV, 
and e' = 42.65A" - 21 (0.925.4 2 - 433 ) MeV, for 
A > 4 (A < 4). The ambient photon distribu- 
tion originates in the thermal emission of the stars 
in the core region of radius R. The average density 
in the region R will reflect both the temperatures 
To and T B due to emission from O and B stars, re- 
spectively, and the dilution resulting from inverse 
square law considerations. Specifically, the photon 
density is 
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where A*o (b) is the number of O (B) stars , i?o(B) 
is the O(B) star average radius, the factor 9/4 
emerges when averaging the inverse square dis- 
tance of an observer from uniformly distributed 
sources in a region R, and 
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Figure 1: Photo-disintegration rates of 56 Fe and 
28 Si, on the Westerlund 2 starlight. 



corresponding to a Bose-Einstein distribution with 
temperature Tb(B) ■ 

In Fig. 1 we show the dependence on the Lorentz 
factor of i?56 and R23 for the stellar ambiance de- 
scribed above. We have taken for the O stars, 
A^o = 12, a surface temperature T = 40000 K, 
and radius Ro — 19 Rq\ for the cooler B stars 
we assign T B = 18000 K, N B = 100, and radius 
Rb = 8 Rq. 

The low energy cutoff on Ra seen in Fig. 1 will 
be mirrored in the resulting photon distribution. 
The ~ E~ 2 energy behavior of the various nuclear 
fluxes will not substantially affect this low energy 
feature. The energy behavior for photons in the 
0.5 — 10 TeV region of the HESS data is a complex 
convolution of the energy distributions of the vari- 
ous nuclei participating in the photo-disintegration, 
with the rate factors appropriate to the eV photon 
density for the various stellar populations. 

Let us define the differential rate dRA/dE' as 
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where da 7 A{^' , E') / dE' is the inclusive differ- 
ential cross section for production of 7-rays from 
disintegration and E' is the energy of the emitted 
photon(s) in the rest frame of the nucleus. As- 
suming the same cosmic ray spectrum as above, 
the emissivity (number/volume/steradian) of 7- 
rays coming from nuclei photo-emission is 



10 



E 

A 



driA 
dE N 



(E N )dE 



N 



dRA dcosOy 

dE„, 

dE^ 7 2 

x <5[£ 7 -7^(1 + cos 7 )] ,(6) 

where E 1 is the energy of the emitted 7-ray in 
the lab and 6 1 is the 7-ray angle with respect to 
the direction of the excited nucleus. Assuming a 
power law with spectral index a for the nuclear 
flux, and approximating the 7-ray spectrum as be- 
ing monochromatic with energy E' A , the emissiv- 
ity becomes [10] 
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where tTa is the mean 7-ray multiplicity for a nu- 
cleus with atomic number A. (Hereafter we take 
Wa = 2). The differential flux at the observer's 
site (assuming there is no absorption) is related to 
to the 7-ray emissivity as 
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where Vdi S is the volume of the source (disintegra- 
tion) region and d is the distance to the observer. 

In Fig. 2 we provide a sample of eyeball fits to the 
gamma ray spectrum, as obtained following a di- 
rect integration of Eq. (7). The fits are for inter- 
esting choices of the spectral index (a) of the nu- 
clei population and for the average energy of the 
photon (in the nuclear rest frame) emitted during 
photo-emission. It is apparent that the A* mech- 
anism can provide reasonable agreement with the 
data, except possibly for the lowest-energy datum. 
This datum may require an alternate mechanism, 
such as electron acceleration or a PION contribu- 
tion [11]. 
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Figure 2: Sample eyeball fits to HESS Wester- 
lund 2 7-ray spectrum for the A* model, with var- 
ious values of the injection 56 Fe nuclei power in- 
dex a and the average energy E' A of emitted de- 
excitation photons. 



The EM and PION pp processes contrast with 
the A* -process in that for them there is either 
no energy threshold (EM) or very small threshold 
0(2771^) (PION pp) in the lab, and so their result- 
ing 7-ray spectra rise monotonically with decreas- 
ing E 1 . In contrast, the A* spectrum is flat with 
decreasing E 1 below about a TeV, as is evident in 
Fig. 2. 

The conversion efficiency from cosmic -rays to 
TeV 7-rays must be quite high in the A* model, as 
with other models. We find that the ratio of power 
in the nuclei flux to that in the O-star winds is 
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where eff is an efficiency which depends on model 
parameters, T ag0 ~ 2 Myr is the age of the star- 
forming region [5], and T CO nt is the containment 
time of the nuclei in the ambient magnetic field. 
For the four parameter sets (a, E'^ A /MeV) listed 
in Fig. 2, we find for the (2.0, 1.5), (2.0, 1.0), 
(1.6, 0.3) and (2.0, 0.3) models, the respective 
efficiency values of eff = 4%, 8%, 20%, and 80%. 
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In summary, we have shown that the observed TeV 
7-rays from Westerlund 2 can be explained by the 
A* model, wherein the TeV 7-rays are the Lorentz- 
boosted MeV 7-rays emitted on the de-excitation 
of daughter nuclei, themselves produced in colli- 
sions of PeV nuclei with a hot ultraviolet photon 
background. There is a specific prediction of a sup- 
pression of the 7-ray spectrum in the region be- 
low 1 TeV - this because of the need to achieve 
Giant Dipole Resonance excitation through colli- 
sion with ~ few eV photons. This suppression is 
not apparent in the lowest-energy Westerlund 2 da- 
tum, but will be probed by the upcoming GLAST 
mission [12]: The flux predicted at - 100 GeV 
from an E~ 2 - 53 extrapolation of the HESS data [2] 
would render the source spectacularly visible in the 
GLAST observation, whereas the A* -model pre- 
dicts a suppression by a factor of more than 2 or- 
ders of magnitude relative to this extrapolated flux. 

We have also calculated the energy requirements 
for the proposed mechanism by integrating over 
the nuclei energy density. We have found that, in 
order to fall within the estimated kinetic energy 
budget of the O stars, the containment time needs 
to be large (of O(10 5 yr)), so that the wind power 
integrates in time to a sufficiently large energy. 
It should be noted that WR 20a may by itself 
contribute several times the wind energy of the 
O stars, thereby lowering our effeciency factor by 
a comparable amount. 
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